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How much matter is there in the universe? Does the universe have the critical 
■ density needed to stop its expansion, or is the universe underweight and destined 
to expand forever? We show that several independent measures, especially those 
utilizing the largest bound systems known - clusters of galaxies - all indicate that 
the mass-density of the universe is insufficient to halt the expansion. A promising 
new method, the evolution of the number density of clusters with time, provides 
the most powerful indication so far that the universe has a sub-critical density. 
We show that different techniques reveal a consistent picture of a lightweight 
universe with only ~ 20-30% of the critical density. Thus, the universe may 

h— » ■ expand forever. 
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Standard models of inflation - how the universe expanded in the beginning - as well as 
general arguments that demand no "fine tuning" of cosmological parameters, predict a fiat 
universe with the critical density needed to just halt its expansion. The critical density, 
1.9 x 10 _29 h 2 g cm -3 (where h refers to Hubble's constant, see below), is equivalent to ~ 10 
protons per cubic meter; this density provides the gravitational pull needed to slow down the 
universal expansion that began with the Big Bang approximately 15 billion years ago, and 
will eventually bring it to a halt. So far, however, only a small fraction of the critical density 
has been detected, even when all the unseen "dark matter" in galaxy halos and clusters of 
galaxies is included. There is no reliable indication so far that most of the matter needed for 
closing the universe does in fact exist. Here we show that several independent observations 
of clusters of galaxies, including the mass-to-light ratio of clusters, the high baryon fraction 
in clusters, and the observed evolution of cluster abundance, all portray a consistent picture 
of a sub-critical universe. 
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1. Weighing Clusters 

Rich clusters of galaxies - families of hundreds of galaxies held together by the gravita- 
tional potential of the cluster - are the most massive bound objects known. Cluster masses 
can be directly and reliably determined using three independent methods: 1) the motion 
(velocity dispersion) of galaxies within clusters reflect the dynamical cluster mass, within a 
given radius, assuming the clusters are in hydrostatic equilibrium (1-3); 2) the temperature 
of the hot intracluster gas, like the galaxy motion, traces the cluster mass (4-6); and 3) 
gravitational lensing distortions of background galaxies can be used to directly measure the 
intervening cluster mass that causes the distortions (7-10). All three independent methods 
yield consistent cluster masses (typically within radii of ~ 1 Mpc ~ 3 x 10 6 light years ), 
indicating that we can reliably determine cluster masses within the observed scatter (~ ± 
30%). 

2. Mass-to-Light Ratio of Clusters 

Let us begin with the simplest argument for a low density universe. The masses of 
rich clusters of galaxies range from ~ 10 14 to 10 15 h _1 M within 1.5h _1 Mpc radius of the 
cluster center (where h = Ho/100 km s _1 Mpc^ 1 denotes Hubble's constant, representing 
the expansion rate of the universe). When normalized by the cluster luminosity, a median 
mass-to-light ratio of M/L B ^ 300±100h in solar units (M /L Q ) is observed for rich clusters, 
independent of the cluster luminosity, velocity dispersion, or other parameters (3,11). (L B 
is the total luminosity of the cluster in the blue band, corrected for internal and Galactic 
absorption.) When integrated over the entire observed luminosity density of the universe, 
this mass-to-light ratio yields a mass density of p m ~ 0.4 x 10 _29 h 2 g cm -3 , or a mass density 
ratio of Q m = p m /p cr it — 0.2 ± 0.07 (where p cr u is the critical density needed to close the 
universe). The inferred density assumes that all galaxies exhibit the same high M/Lb ratio 
as clusters, and that mass follows light on large scales. Thus, even if all galaxies have as 
much mass per unit luminosity as do massive clusters, the total mass of the universe is only 
~ 20% of the critical density. If one insists on esthetic grounds that the universe has a 
critical density (Q m = 1), then most of the mass of the universe has to be unassociated with 
galaxies (i.e., with light). On large scales (^ 1.5 h" 1 Mpc) the mass has to reside in "voids" 
where there is no light. This would imply, for Q m = 1, a large bias in the distribution of 
mass versus light, with mass distributed considerably more diffusely than light. 

Is there a strong bias in the universe, with most of the dark matter residing on large scales, 
well beyond galaxies and clusters? A recent analysis of the mass-to-light ratio of galaxies, 
groups, and clusters by Bahcall, Lubin and Dorman (11) suggests that there is not a large 
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bias. The study shows that the M/L B ratio of galaxies increases with scale up to radii of 
R ~ 0.2 fr 1 Mpc, due to very large dark halos around galaxies (see also 12,13). The M/L 
ratio, however, appears to flatten and remain approximately constant for groups and rich 
clusters from scales of ~ 0.2 to at least 1.5 h _1 Mpc and possibly even beyond (Fig.l). The 
flattening occurs at M/Lb — 200 — 300h, corresponding to Vt m ~ 0.2. (An M/Lb ~ 1350h 
is needed for a critical density universe, Q m = 1.) This observation contradicts the classical 
belief that the relative amount of dark matter increases continuously with scale, possibly 
reaching Vl m = 1 on large scales. The available data suggest that most of the dark matter 
may be associated with very large dark halos of galaxies and that clusters do not contain a 
substantial amount of additional dark matter, other than that associated with (or torn-off 
from) the galaxy halos, plus the hot intracluster gas. This flattening of M/L with scale, 
if confirmed by further larger-scale observations, suggests that the relative amount of dark 
matter does not increase significantly with scale above ~ 0.2 h -1 Mpc. In that case, the 
mass density of the universe is low, Vt m ~ 0.2 — 0.3, with no significant bias (i.e., mass 
approximately following light on large scales). 



3. Baryons in Clusters 

Clusters contain many baryons, observed as gas and stars. Within 1.5h _1 Mpc of a rich 
cluster, the X-ray emitting gas contributes ~ 6h -L5 % of the cluster virial mass (14-16). 
Stars contribute another ^ 4%. The baryon fraction observed in clusters is thus: 

n b /n m £ o.o6h" 15 + o.o4 (i) 

Standard Big Bang nucleosynthesis limits the baryon density of the universe to (17-18): 

n b ~ 0.017h" 2 (2) 

These facts suggest that the baryon fraction observed in rich clusters (eq.l) exceeds that of 
an Q m = 1 universe (f2&/(f2 m = 1) ~ 0.017h -2 ; eq. 2) by a factor of ^ 3 (for h ^ 0.5). 
Since detailed hydrodynamic simulations (14,16) show that baryons do not segregate into 
rich clusters, the above results imply that either the mean density of the universe is lower 
than the critical density by a factor of ^ 3, or that the baryon density is much larger than 
predicted by nucleosynthesis. The observed high baryonic mass fraction in clusters (eq.l), 
combined with the nucleosynthesis limit (eq.2), suggest (for h ~ 0.5 — 1) 



Q m ~ 0.2 ±0.1. 



(3) 
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4. Evolution of Cluster Abundance 

In a recent study by Bahcall, Fan and Cen (19,20) we show that the evolution of the num- 
ber density of clusters as a function of cosmic time (or redshift) provides a powerful constraint 
on fl m (19-22). The growth of high-mass clusters from initial Gaussian fluctuations depends 
strongly on the cosmological parameters Q m and <x 8 (where eg is the root-mean-square mass 
fluctuation on 8 h _1 Mpc scale; 23-27). In low-density models, density fluctuations evolve 
and freeze out at early times, thus producing only relatively little evolution at recent times 
(z < 1). In an Q m = 1 universe, the fluctuations start growing more recently thereby pro- 
ducing strong evolution in recent times; a large increase in the abundance of massive clusters 
is expected from z ~ 1 to z ~ 0. The evolution is so strong in Q m = 1 models that finding 
even a few Coma-like clusters at z > 0.5 over ~ 10 3 deg 2 of sky contradicts an Q m = 1 
model where only ~ 10~ 2 such clusters would be expected (when normalized to the observed 
present-day cluster abundance). The evolution of the number density of Coma-like clusters 
was recently determined from observations and compared with cosmological simulations (19- 
21). The data show only a slow evolution of the cluster abundance to z ~ 0.5, with ~ 10 2 
times more clusters observed at these redshifts than expected for Vl m = 1. The results yield 

a rn ~ 0.3 ±o.i. 

The evolutionary effects increase with cluster mass and with redshift. The existence 
of the three most massive clusters observed so far at z ~ 0.5 — 0.9 places the strongest 
constraint yet on Q m and a 8 . These clusters (MS0016+016 at z = 0.55, MS0451+03 at 
z = 0.54, and MS1054-03 at z = 0.83, from the Extended Medium Sensitivity Survey, EMSS, 
28) are nearly twice as massive as the Coma cluster, and have reliably measured masses 
(including gravitational lensing masses, temperatures, and velocity dispersions; 3,9,29-32). 
These clusters posses the highest masses (^,8 x 10 14 h _1 M within 1.5 h _1 comoving Mpc 
radius), the highest velocity dispersions (^ 1200 km s _1 ), and the highest temperatures (^ 8 
kev) in the z > 0.5 EMSS survey. The existence of these three massive distant clusters, even 
just the existence of the single observed cluster at z — 0.83, rules out Gaussian Vl m = 1 
models for which only ~ 10~ 5 z ~ 0.8 clusters are expected instead of the 1 cluster observed 
(or ~ 10~ 3 z > 0.5 clusters expected instead of the 3 observed). (See Bahcall & Fan (29)). 

In Figure 2 we compare the observed versus expected evolution of the number density of 
such massive clusters. The expected evolution is based on the Press-Schechter (23) formal- 
ism that describes the growth of structure in a hierarchical universe with standard initial 
Gaussian density fluctuations; this formalism agrees well with direct numerical cosmological 
simulations (20,26). The expected evolution is shown for different Q m values (each with 
the appropriate normalization <7g that satisfies the observed present-day cluster abundance, 
a 8 ~ 0.5f2~ a5 ; 26,33). The model curves range from Q m = 0.1 (<x 8 — 1-7) at the top of the 
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figure (flattest, nearly no evolution) to Vt m — 1 (<r 8 ~ 0.5) at the bottom (steepest, strongest 
evolution). The difference between high and low Q m models is dramatic for these high mass 
clusters: Q m = 1 models predict ~ 10 5 times less clusters at z ~ 0.8 than do fl m ~ 0.2 mod- 
els. The large magnitude of the effect is due to the fact that these are very massive clusters, 
on the exponential tail of the cluster mass function; they are rare events and the evolution 
of their number density depends exponentially on their "rarity", i.e., depends exponentially 
on a 8 2 oc Q m (20,23,29). The number of clusters observed at z ~ 0.8 is consistent with 
Q m ~ 0.2, and is highly inconsistent with the ~ 10~ 5 clusters expected if Vl m = 1. The 
data exhibits only a slow, relatively flat evolution; this is expected only in low -fl m models. 
tt m = 1 models have a ~ 10~ 5 probability of producing the one observed cluster at z ~ 0.8, 
and, independently, a ~ 10~ 6 probability of producing the two observed clusters at z ~ 0.55. 
These results rule out fl m = 1 Gaussian models at a very high confidence level. The results 
are similar for models with or without a cosmo logical constant. The data provide powerful 
constraints on Q m and a 8 : Q m = 0.2t° f and a 8 = 1.2 ±0.3 (68% confidence level) (29). The 
high <7g value for the mean mass fluctuations indicates a nearly unbiased universe, with mass 
approximately tracing light on large scales (since the galaxy fluctuations, which represent 
the light, exhibit a similar value of a 8 (galaxy) ~ 1). This conclusion is consistent with the 
suggested flattening of the observed M/L ratio on large scales (Fig. 1). 

In Figure 3 we summarize the four independent Vt m (cr 8 ) constraints obtained from the 
cluster results discussed above: 1) the present-day cluster abundance constraint (26, 33) 
Vt ^ ~ 0.5/(T 8 ; 2) the high-redshift (z ~ 0.5 — 0.9) cluster abundance constraint (29); [the 
overlap of the z ~ and z ~ 0.5 — 0.9 abundance constraints of (1) and (2) yields the cluster 
evolution constraint discussed above]; 3) the fl m derived from the high baryon fraction in 
clusters; and 4) the Q m obtained from cluster masses. The results are all consistent with 
each other for Q m = 0.2 ± 0.1 and a 8 = 1.2 ± 0.2 (ler level). Q m = 1 models are highly 
incompatible with these results (^ 10~ 6 probability). 



5. Summary 

We have shown that several independent observations of clusters of galaxies all indicate 
that the mass-density of the universe is sub-critical: Q m ~ 0.2 ± 0.1. A summary of the 
results, presented in Fig. 3, is highlighted below. 

1. The mass-to-light ratio of clusters of galaxies and the suggested flattening of the 
mass-to-light ratio on large scales suggest Vl m ~ 0.2 ± 0.1. 

2. The high baryon fraction observed in clusters of galaxies suggests Vl m ~ 0.2 ± 0.1. 
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3. The weak evolution of the observed cluster abundance to z ~ 1 provides a robust 
estimate of Q m ~ 0.2l[J;J 5 , valid for any Gaussian models. An Q m = 1 Gaussian universe is 
ruled out as a ^ 1CT 6 probability by the cluster evolution results (Fig. 2-3). 

4. All the above-described independent measures are consistent with each other and indi- 
cate a low-density universe with Q m ~ 0.2 ±0.1 (Fig. 3). VL m = 1 models are ruled out by the 
data. While non-Gaussian initial fluctuations, if they exist, will affect the cluster evolution 
results, they will not affect arguments (1) and (2) above. Gaussian low-density models (with 
or without a cosmological constant) can consistently explain all the independent observations 
presented here. These independent cluster observations indicate that we live in a lightweight 
universe with only ~ 20% — 30% of the critical density, Thus, the universe may expand 
forever. 

We thank D. Eisenstein, J. P. Ostriker, P. J. E. Peebles, D. N. Schramm, and D. N. Spergel 
for helpful discussions. This work was supported in part by NSF grant AST93-15368. 
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Figure l.The dependence of mass-to-light ratio, M/Lb, on scale, R, for average spiral galaxies 
(blue symbols), elliptical galaxies (green), and groups and clusters (red). (From Bahcall, 
Lubin and Dorman 1995)(11). The large scale point at ~ 15h _1 Mpc represents Virgo 
cluster infall motion results (11). The location of Q m = 1 and Q m = 0.3 are indicated by the 
horizontal lines. A flattening of M/Lb is suggested at fl m ~ 0.2 ± 0.1. 
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Figure 2. Evolution of the number density of massive clusters as a function of redshift: 
observed versus expected (for clusters with mass ^8x 10 14 h _1 M Q within a comoving radius 
of 1.5h _1 Mpc). From Bahcall and Fan 1998 (29). The expected evolution is presented for 
different VL m values by the different curves. The observational data points (see text) show 
only a slow evolution in the cluster abundance, consistent with fl m ~ 0.21q^ 5 . Models with 
Vl m = 1 predict ~ 10 5 fewer clusters than observed at z ~ 0.8, and ~ 10 3 fewer clusters than 
observed at z ~ 0.6. 
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Figure 3. Constraining the mass-density parameter, Q m , and the mass fluctuations on 8h _1 
Mpc scale, as, from several independent observations of clusters: cluster dynamics (blue 
band); baryon fraction in clusters (pink); present-day cluster abundance (z ~ 0; green); and 
cluster abundance at redshift z ~ 0.7 (yellow). (The latter two abundances yield the cluster 
evolution constraints shown in Fig. 2; see text). All these model-independent observations 
converge at the allowed range of Q m = 0.2±0.1 and as = 1.2±0.2 ( 68% confidence level; red). 
The dotted lines illustrate the mean microwave fluctuations constraints, based on the COBE 
satellite results, for a Cold-Dark-Matter model with h = 0.7 (with and without a cosmological 
constant, denoted as LCDM and OCDM respectively. Both models are consistent, within 
their uncertainties, with the best-fit Q m — a 8 regime of the cluster observations). 



